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Abstract 

The dynamics of fluctuating radially growing interfaces is approached using the formalism of 
stochastic growth equations on growing domains. This framework reveals a number of dynamic 
features arising during surface growth. For fast growth, dilution, which spatially reorders the 
incoming matter, is responsible for the transmission of correlations. Its effects include the erasing 
of memory with respect to the initial condition, a partial regularization of geometrically originated 
instabilities, and the restoring of universality in some special cases in which the critical exponents 
depend on the parameters of the equation of motion. All together lie in the basis of the preservation 
of the Family- Vicsek scaling in radial interfaces, which is thus a direct consequence of dilution. 
This fast growth regime is also characterized by the spatial decorrelation of the interface, which 
in the case of radially growing interfaces naturally originates rapid roughening and polyfractality, 
and suggests the advent of a self-similar fractal dimension. The center of mass fluctuations of 
growing clusters are also studied, and our analysis supports a strong violation of the Family- Vicsek 
scaling by the surface fluctuations of the Eden model. Consequently, this model would belong to 
a dilution-free universality class. 

PACS numbers: 64.60.Ht, 05.10.Gg, 05.40.-a, 68.35. Ct 
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I. INTRODUCTION 



The study of fluctuating interfaces has occupied an important place within statistical 
mechanics in recent and not so recent times. The origins of this interest are practical, due 
to the vast range of potential applications that this theory may have, and theoretical, as 
some of the universality classes discovered within this framework are claimed to play an 
important role in other areas of physics [1]. While the great majority of works on this topic 
has concentrated on strip or slab geometries, it is true that at the very beginning of the 
theoretical studies on nonequilibrium growth one finds the seminal works by Eden, focused 
on radial shapes To a certain extent, the motivation of considering radial forms is 

related to biological growth, as for instance the Eden model can be thought of as a simplified 
description of a developing cell colony. The Eden and other related discrete models have 
been computationally analyzed along the years, and the results obtained has been put in 
the context of stochastic growth theory, see for instance {4] and references therein. 

Apart from the interest in modelling, there is a genuine theoretical motivation in un- 
derstanding the dynamics of growing radial clusters. The Eden model is actually a sort of 
first passage percolation [sl, and the scahng limit of percolation models has been studied 
by means of field-theoretic approaches [o] and stochastic processes like Schramm-Loewner 
evolution [7|]. A natural theoretical question to be answered is in which cases the Family- 
Vicsek scaling Sj], basic to describe planar growth processes, is able to capture the behavior 
of the surface fluctuations of growing radial clusters. 

The use of stochastic differential equations, very much spread in the modelling of planar 



growth profiles, has been not so commonly emp 



works constitute an exception to this rule 



oyed in the case of radial growth. A series of 
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15|, as they proposed a partial 



differential equation with stochastic terms as a benchmark for analyzing the dynamics of 
radial interfaces. Because studying this sort of equations is complicated by the nonlinearities 
implied by reparametrization invariance, a simplified version in which only the substrate 
growth was considered was introduced in [l^. Already in this case it was apparent that 
for rapidly growing interfaces dilution, which is responsible for matter redistribution as 
the substrate grows [l3], propagates the correlations when large spatiotemporal scales are 
considered. It is also capable of erasing the memory effects that would otherwise arise, let 
us show how. In [l^ we considered the linear equation for stochastic growth on a growing 
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domain 



dth = -D M'^h -^h+ ^Ft^-' + '''' ax, 



(1) 



where the domain grows following the power law if:'*', 7 > is the growth index and —{d'y/t)h 
is the term taking into account dilution IGj. Its Fourier transformed version, for > 1, is 



dhn 
~dt 



-D 



d'-y , 



Lq t 



(2) 



This equation can be readily solved for 7 > 1/^ and in the long time limit 



hn{t) = (t/to)-"^ 



exp 



Dto TT^Inl"^' 



c 



hn{to) + {t/h)-"^ 



to 



to J 



dj/2 



Ur)dT, (3) 



and so the dependence on the initial condition tends to zero as a power law for long times. 
This is, as mentioned, one of the consequences of dilution. If we considered the dilation 
transformation x —>■ {t/tQ)''x we would find again Eq. ([T]) but this time without the dilution 
term. The solution now becomes 



h„(t) = exp 



Dto vr'^kl'^' 



hn{to) + 



(4) 
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and so the dependence on the initial condition remains for all times. In the first case the 
long time solution becomes spatially uncorrelated, and in the second one only part of the 
initial correlations survive. As an abuse of language, we will talk about decorrelation in 
both cases. The memory effects that affect the solution in the no-dilution (or dilation) 

nn 

situation separate its behavior from the one dictated by the Family- Vicsek ansatz |16l . |18| . 
For 7 < the memory effects and the corresponding dependence on the initial condition 
disappear exponentially fast for long times as a consequence of the effect of diffusion. 

Dilution is also the mechanism that controls the amount of matter on the interface. Pure 
diffusion on a growing domain is described by the equation 



to 



t 
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d'j . 



d,h = D{ - ] V'h -h 



t 



(5) 



in Eulerian coordinates x E [0, Lq] (see |16|) and where dilution has been taken into account. 
The total mass on the surface is conserved 



h{y,t)dy 



d'j 



Lo 



h{x, t)dx 



Lo 



h{x,to)dx, (6) 
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where y = [L(t)/Lo]x denotes the Lagrangian coordinates. In the no-dilution situation we 
find 

■ j h{y,t)dy=i—j J ■J h{x,t)dx = i—j J ■ j h{x,to)dx. 

(7) 

This second case is pure dilation, which implies that not only the space grows, but also the 
interfacial matter grows at the same rate, in such a way that the average density remains 
constant. Note that this process of matter dilation, as well as the spatial growth, are 
deterministic processes. These calculations show that both dilution and dilation dynamics 
are physically motivated and have a number of measurable differences. It is worth remarking 
here that all previous works except [l^ have exclusively considered dilation dynamics. Even 
in the different field of react ion- diffusion dynamics in which the dilution term was derived, 



the focus was on the limit in which it was irrelevant |l7l |. 

This work is devoted to further explore the consequences of dilution, dilation and decorre- 
lation, and their effects on scaling of radial interfaces. We will use in cases radial stochastic 
growth equations, which may show up instabilities , and explore the interplay of dilution 
with them. In other cases, when instabilities do not play a determinant role and for the sake 
of simplicity, we will consider stochastic growth equations on growing domains. 

II. RADIAL RANDOM DEPOSITION 

In order to construct radial growth equations one may invoke the reparametrization 
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invariance principle [19|,[2^, as has already been done a number of timesjo],!! 
In case of white and Gaussian fluctuations, the d— dimensional spherical noise is given by 

1 



=aO,t), ^e(^,t)) = 0, (8) 

^{e,me,t)) = e5{e-e')5{t-t'), (9) 



''9 



where g = det{gij) = det(9jr ■ dj-r) is the determinant of the metric tensor. Under the 
small gradient assumption |Vg r| ^ r one flnds g ~ J'{r,6y, where J' is the Jacobian 
determinant of the change of variables from the Cartesian representation {x, h) to the polar 
representation {6,r). We also have the factorization J{r^9Y = r'^'^J{6)'^, where J is the 
Jacobian evaluated at r = 1. 
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The simplest growth process is possibly the radial random deposition model. If the growth 
rate is explicitly time dependent, then the growth equation reads 

in the absence of dilution. Here r{9, t) is the radius value at the angular position 6 and time 
t, F > is the growth rate, 7 > is the growth index, d is the spatial dimension and ^ is a 
zero mean Gaussian noise, which correlation is given by 



^^{e,t)ae',s)) = e6{e~e')6{t-s). (ii) 

The equation for the first moment can be easily obtained 

dt (r) = F-ft'r-\ (12) 
due to the Ito interpretation, and integrate it to get 

[r{e,t)) = Ft\ (13) 



where we have assumed the radially symmetric initial condition r{6,tQ) = Ft^ and to ^ 
t is the absolute origin of time. It is difficult to obtain more information from the full 
equation (|TOl) . so we will perform a perturbative expansion. We assume the solution form 

r{e,t) = R{t) + yrep^{0,t), (14) 

I — I 

where the noise intensity e will be used as the small parameter [21]. Substituting this solution 
form into Eq. flTU]) we obtain the equations 

dtR = F^t^-\ (15) 

where ^ = i/er]. These equations have been derived assuming ^/e <^ Ff^, a condition much 
more favorable (the better the larger 7 is) than the usual time independent ones supporting 



small noise expansions 



2l|. The solution to these equations can be readily computed 

R{e,t) = Ft\ (17) 
Pi{0,t)) = 0, (18) 



p^ie,t)p^{e\s)) = imm{t, s})'-^' - t]-^'] (19) 



J{9) 



if 7*^ 7^ 1 and where we have assumed a zero value for the initial perturbation. If 7(i = 1 the 
correlation becomes 



p,{e\t)p,{e\s)) = ^\xi 



min{t, s} 



(20) 

m 



Here i? is a deterministic function and pi is a zero mean Gaussian stochastic process that 
is completely determined by the correlations given above. The long time behavior of the 
correlations, given by the condition t, s ^ to, is specified by the following two-times and 
one-time correlation functions 

pr(e,t)p,(e',s)) = -^(mmli.ii})'-^''^*^^, (21) 



if 7^ > 1, 



p,mp,{e\s)) = i^ln(min{t,.})^i^^, (23) 

/ F j{e) 

pr{e,t)p,{e',t)) = ^ln(t)^^^;^, (24) 



if 7(i = 1, and finally 



when '-yd > 1. In this last case the correlation vanishes in the limit to oo. Note that the 
reparametrization invariance principle is not able to capture dilution effects and it reproduces 
pure dilation dynamics. 

In order to introduce dilution in the radial case we may use the following functional 
definition which transforms Eq. ffTOj) into 



Otr = F'yV r + —Tj^ — , (26) 

whose first moment can be exactly calculated again taking advantage of the Ito interpretation 
of the noise term, yielding 
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Performing as in the former case the small noise expansion r = R + \f^P\ we find again 
R= (r). The perturbation obeys the equation 

and so the perturbation has zero mean and its long time correlation is given by 

l^p,{Q^)p,{Q',s)) = ^i^^V/l) "^M^,Omax{.,t}-^'^^^^^, (29) 

a result that holds uniformly in 7. Note that the structure of the temporal correlation 
is different when the effect of dilution is considered and when it is not for all 7 > 0. 
For instance, the characteristic length scale corresponding to a given angular difference is 
A = max{s, t}'^\6 — 6'\ when dilution is present, and A = min{s, t}'^\6 — Q'\ in the absence of 
dilution. One already sees in this example that the lack of dilution causes the appearance 
of memory effects on the growth dynamics. The first order correction in the small noise 
expansion p\ is always a Gaussian stochastic process; an attempt to go beyond Gaussianity 
by deriving the second order correction is reported in appendix El 



III. RANDOM DEPOSITION AND DIFFUSION 



Our next step, in order to approach more complex and realistic growth processes, is 
to add diffusion to a random deposition equation of growth. This sort of equations may 
be derived using reparametrization invariance as in [131. Following this reference and the 
former section, we perform a small noise expansion and concentrate on the equation for the 
Gaussian perturbation. In this section we will consider a number of cases which do not show 
instabilities, and the study of these will be postponed to the next one. The equation for the 
perturbation in ci = 1 is 15 1 

where Ag is a fractional differential operator of order and dilution has not been considered. 
The dynamics for ( > d, which in turn implies that in the linear case the growth exponent 
(3 > and the interface is consequently rough, has been already considered in [l^; herein we 
move to studying the marginal case ( = d, which turns out to have interesting properties. 
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The case ( < d is not so interesting as it corresponds to flat interfaces; an analogous 
calculation to the corresponding one in for 7 = 1 and C < 1 shows 



when 



s ^ 00, 



(31) 



independently of the value of to- 
If C = 7 = 1 the correlation reads 

1 



{p{e,t)p{e',s)) 



A-nD 



In 



{ts) 



D/F 



{s/t)D/P + {tls)DIP - 2 cos(^ - 6' 



The one time correlation adopts the form 



{p{e,t)p{e\t)) 



1 



47rL> 



-In 



2 - 2 COS - 9') 



that reduces to 



{p{e,t)p{e\t)) 



-In 



t 



(32) 



(33) 



(34) 



2tiF"' \ \0-e'f''^ 

when we consider local in space dynamics, this is, in the limit 9 9'. Note that this 
result allows us to deflne the local dynamic exponent zioc = F/D E (0, 00), which depends 
continuously on the equation parameters F and D, and is thus nonuniversal, as we noted 
in In terms of the arc-length variable i — i' = t{9 — 9') we flnd 



{p{^,t)p{^\t)) 



D 



-1 



27r 



-In 



t 



■nF/{D+F) 



(35) 



where the dynamical exponent in terms of the arc-length variable ze = F/{D + F) G (0, 1) 
is again nonuniversal. If we take into account dilution Eq. fl30|) transforms to 

dtp = §-Agp - + -^V{0, t). (36) 



(37) 



The solution has zero mean and its correlation is given by 



47ri' 

(min{s, t} I max{s, t)Y"^^l^ 



2-n{F + D) 



3? <^ d^'-''\F^ 



' D' £)' \max{s,t} 



where 3fJ(-) denotes the real part and 2-Fi(-, ■; ■; ■) is Gauss hypergeometric function [22|]. This 
correlation, for s = t and for small angular scales 9 ^ 9\ becomes at leading order 

-1 



{p{9,t)p{9\t)) 



27tD 



In(l^-^'l) 



(38) 



which is time independent, and for the arc-length variable 



(p(£,t)p(f,t)) 



1 



-In 



t 



(39) 



2ttD"' 

for which the planar scaling and the universal dynamical exponent z = 1 are recovered, see 
Eq. (C5) in [15]. This is yet another example, this time of a different nature, of how dilution 
is able to restore the Family- Vicsek ansatz [igI, [3] ■ 

If ( = 1 and 7 < 1 we find the following correlation function 



[minjif:, s}] 



1-7 



27rF(l-7) AnD 
—2 exp 

When t = s we get 

{p{e,t)p{e',t)) 



In 1 + exp 
D 



2D 



F(l-7) 



a-7| 



,1-7 
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.1-71 



1 



27rF(l - 7) AttD 
and considering local spatial dynamics we arrive at 

t^-^ 1 



In [2-2 cos (e - e')] 



{p{e,t)p{e',t)) 



-ln(l^^-^^'l) 



1 



-In 



+1 — t 



\0 _ 5)/|F(l-7)/L> 



(40) 



(41) 



(42) 



27rF(l-7) 2nD " " 27rF(l - 7) 

expression that does not allow to define a local dynamic exponent, or alternatively zioc = 
due to the exponentially fast spreading of the correlations. These last three expressions 
contain two clearly different terms. The first one is the zeroth mode component of the 
correlation, which does not achieve long time saturation. The second term is the nontrivial 
stationary part of the correlation generated along the evolution. As can be seen, both spatial 
and temporal correlations are generated. 

When the dilution term is taken into account we find the correlation 



{p{e,t)p{e',s)) 



min{t, s}[max{t, s}] 



27rF(7 + 1) AnD 
—2 exp 

When t = s we get 

{p{e,t)p{e',t)) 



In 1 + exp 
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27rF(7+l) AttD 
and considering local spatial dynamics we arrive at 
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(45) 
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and we see that as in the former case, both prefactor and exponent are modified, but the still 
exponentially fast propagation of correlations implies an effective local dynamical exponent 
zioc = 0. Note that for 7 > 1 a radial random deposition behavior for large spatial scales is 
recovered. 

Now we move onto the two-dimensional setting. As in the one-dimensional case we focus 
on the marginal situation d = ( = 2, which leads us to denominate this sort of equations 
as spherical Edwards- Wilkinson (EW) equations, and < 7 < 1/2, as greater values of the 
growth index lead again to decorrelation. The straightforward generalization of Eq. (130|) is 



dtp 



K 



1 



where the noise is a Gaussian random variable of zero mean and correlation given by 

{^e, 0, me', 0', s)) = 6{e - e')6{<p - <p')6{t - s). 



(46) 



(47) 



In this case, if 7 < 1/2, the random variable p is a zero mean Gaussian process whose 
correlation is given by 



00 / 

EE 
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(=1 m=—l 
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exp 



47rF2(l - 27) 



r!„(e,0)il(^',0') 



(4J 



F2(l - 27) I 

where the expansion has been performed on the spherical harmonics basis Y^{9,(f)). If 
7 = 1/2 then p becomes a zero mean Gaussian random variable with the new correlation 

In [min(t, s)] 



°° ^ ( -\\m 

^ ^ 2K{1 + P 

1=1 m=-l ^ 



{p{e,Mp{e'A',s)) 

K(l+l^)/F'^ 



min(s, t) 
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yl 



„(^,0)ll(^',,-' 



)• 



(49) 



It is clear that these correlations are again composed of two different terms, the first one 
associated with the / = mode never saturates, and the second one associated with the rest 
of modes I > 0, which saturates and is responsible of a non-trivial spatial structure. 
Taking into account dilution we find for 7 < 1/2 the correlation 

/ //I , N ,/ min(t, s) [maxft, s)!"^''' 



47rF2(27 + 1^ 



00 I 

EE 

(=1 m=—l 



2K{1 + /2) 



exp 
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F2(l - 27) 
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(50) 
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and for 7 = 1/2 

oo ( / ^ 



(p(^,0,t)p(^',0',s)) 



mm{s,tj 



max(s, t) 



2F2 + 2K(P + /) 

(51) 

In the two dimensional situation we see that dilution also has a measurable effect, which is 
more pronounced in the critical 7 = 1/2 case. For this value all the modes in the correlation 
saturate and contribute to create a stationary spatial structure, as in the one-dimensional 
setting. It is difficult to establish more comparisons among both dimensionalities, as the 
infinite sums that were explicit in ci = 1 become much more involved in d = 2, due to the 
double series containing the spherical harmonics. We however conjecture that the modifi- 
cation of the scaling properties due to effect of dilution in two dimensions is similar to the 
one explicitly observed in one dimension. 

IV. INSTABILITIES 

A spherical EW equation derived from the geometric principle of surface minimization 



was introduced in 

dtr = K 



15|. The corresponding equation for the radius r{9,(p,t) reads 



der ^ 9|r ^ _ 2 



+ F^r-' + ^==ae,<P,t). (52) 



Performing the small noise expansion r{6,(f),t) = Ft + p{6,(j),t) we find a linear equation 
which differs from Eq. (H6l) in that it has a destabilizing term coming from the fourth term 
in the drift of Eq. flS21) . see [15]. In this reference one can see that in the absence of dilution 
the I = mode is unstable and the / = 1 modes are marginal while the rest of modes is 
stable. The effect of this sort of geometrically originated instability on the mean value of the 
stochastic perturbation and alternative geometric variational approaches that avoid it can 
be seen in herein we will concentrate on its effect on correlations. Its effect on mean 
values can be easily deduced from them. 

In the long time limit and provided 7 < 1/2, the perturbation is a Gaussian process 
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whose correlation is given by 

{p{e,4>.t)p{e\<p\s)) 
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47rF2(l - 27) 



cos(6') cos(6'') + cos(^ - (j)') sin(6') sin(6'')] + 



00 I 

1=2 m=-l 



-11 



2K{P + 1-2) 



exp 



K{1^ + 1-2) I , 



F2(l - 27) 
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(53) 



If 7 = 1/2 the correlation shifts to 
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(54) 



In these cases the modes characterized by Z = and Z = 1 do not saturate, and the rest of 
the modes / > 1 saturate and create a non-trivial spatial structure. When 7 < 1/2 the I = 1 
modes grow in time as a power law with the exponent 1 — 27, while the / = mode grows 
exponentially fast. When 7 = 1/2 the / = 1 modes grow logarithmically and the / = mode 
grows as a power law with the non- universal exponent 4K/F'^. 

When we consider the effect of dilution, and for 7 < 1/2, we find the correlation 



{pi9,cj>,t)pi9',cj>',s)) 
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For 7=1/2 the correlation reads 



F2(l - 27) 



^x-27_ ^1-271 
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where 



(2F2 -4ir)-i (min{s,t}/max{s,t})^"^^/-^' if F^>2K, 
ln{mm{s, t})/F^ if F^^2K, (57) 



(AK - 2F2)-i (ts/tl) 



if 



< 2K, 
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where to is the absolute origin of time. 

Contrary to what happens in the stable case, Eq. P6|) . in the unstable case with no 
dilution, Eq. (13^ . the / = mode is unstable, showing an exponential growth, and the / = 1 
modes shows an algebraic increase with the universal exponent 1 — 27, provided 7 < 1/2; 
the rest of modes is stable. The marginal value of the growth index 7 = 1/2 translates into 
a power law increase of the / = mode with a non- universal exponent, while the / = 1 
modes grow logarithmically; the rest of modes is again stable. It is clear that dilution has 
a stabilizing effect. Indeed, for 7 < 1/2 the / = mode is unchanged, but the / = 1 modes, 
which still grow in time, experience a lost of memory effects. In the critical 7 = 1/2 situation 
the dilution effects are stronger. The / = 1 modes, which formerly grew logarithmically, now 
become stable; the / = mode, which formerly showed an algebraic growth, now shows 
(non-universal) algebraic or logarithmic grow, or even saturation, depending on the relation 
among the values of the parameters of the spherical EW equation. In any case, even that 
of algebraic growth, this growth is always slower than in the no dilution situation. Stable 
modes saturate contributing to create a non-trivial spatial structure in the whole range 
7 < 1/2. 

In summary, the effect of dilution is weakly stabilizing in the subcritical case, while 
stronger and more identifiable in criticality. Of course, the supercritical situation is charac- 
terized by an effective random deposition behavior in the large spatial scale. 

V. INTRINSICALLY SPHERICAL GROWTH AND RAPID ROUGHENING 

It is necessary to clarify the role of the diffusivity index (. We have defined it as the order 
of the fractional differential operator taking mass diffusion into account, and so far we have 
referred to it as the key element triggering decorrelation. This has been an abuse of language 
because we have assumed that the negative power of the radius (or its mean field analog 
Ff^ - what really matters is the resulting power of the temporal variable) preceding this 
differential operator was exactly (. This would not be the case if the diffusion constant were 
time or radius dependent, but also in some other cases as the Intrinsically Spherical (IS) 
equation derived from geometric variational principles in [15]. This equation was obtained 
as a gradient flow pursuing the minimization of the interface mean curvature, and then 
linearizing with respect to the different derivatives of the radius as given by the small gradient 
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assumption [15|. It is termed "intrinsically spherical" because it has no planar counterpart, 
as the nonlinearity becomes fundamental in any attempt to derive such a gradient flow in 
the Cartesian framework 23|. Note the similarity of this with other equations classical in 



this context, as the EW equation is a gradient flow which minimizes the surface area and 
the MuUins-Herring equation minimizes the interface square mean curvature the IS 
equation, as mentioned, minimizes the interface mean curvature. It reads 



i5i 



dfT = K 



+ 



dlr 



+ 



der 



1 



(58) 



r3sin2(^) ' r3tan(0) r2j ' ^ ' r^^J^ 

and so C = 2 in this case, but however one finds a factor in front of the diffusive differential 
operator, instead of the r"^ factor characteristic of the EW equation. This difference will 
have a number of measurable consequences, as we will show in the following. The equation 
for the stochastic perturbation reads in this case 



(59) 



which reveals that the critical value of the growth index 7 = 1/3; a faster growth leads 
to decorrelation. This is the first but not the unique difference with respect to the EW 
equation. To find out more we will first put things in a broader context. 

A more general equation for radial growth, after introducing dilution, would be 



dtr = -4|V|^r - — r + F-ft^'^ + 
r° t 



(60) 



which defines the damping index 5, differing from the diffusivity index ( in general; note that 
Eq. ( I6OI) has left aside the instability properties of the IS equation, which are analogous to 
those of the EW equation, and would add nothing to last section discussion. For simplicity 
we will focus on values of the damping index fulfilling S > (. This equation can be treated 
perturbatively for small e following the previous sections procedure and by introducing the 
hyperspherical harmonics Yf^{6), which obey the eigenvalue equation 24 1 



v'Yi^ie) = -i{i + d- i)Yi"'{e), 



(61) 



where the vector m represents the set of {d — 1) indices. The fractional operator acts on the 
hyperspherical harmonics in the following fashion 



\V\'^Yr{9) = + d- 1)]'^/'1^™(^). 



(62) 
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The hyperspherical noise is Gaussian, has zero mean and its correlation is given by 



7]{e,t)r]{e',t')) = 6{e-e')6{t-t'). 



(63) 



It can be expanded in terms of hyperspherical harmonics 



J{9) l,rn 



(64) 



and the amplitudes are given by 



vfit) = [ vio,t)Yrie)^)d9, 



(65) 

and so they are zero mean Gaussian noises whose correlation is given by 

[vfitW'it') )=5it- t')5l,,6rn,rn', (66) 

where the overbar denotes complex conjugation. Note that the amplitudes are in general 
complex valued. They obey the linear equation 



dpi 



K 



■[/(/ + rf-l)pp[ 



C/2 m 7*^ r» 



1 



(67) 



From this equation it is clear that the critical value of the growth index is 7 = 1/6, and a 
faster growth leads to decorrelation. 



It is convenient to move to a growing hypercubic geometry as in 
different quantities 



16j in order to calculate 



x,t), 



(68) 



since this change simplifies calculations without modifying the leading results. Our goal is 
finding the growth and auto-correlation exponents, as this last one is a good quantity to 



measure decorrelation 



16(1 . In order to calculate the temporal correlations we need to con- 



sider the short time limit, where the growth exponent (3 becomes apparent. The propagator 
of Eq.§Ei) is 



Gnit) = { f 
^0 



exp 



n 



'^TT^D tft^--^^ - to 



4 



1-7(5 

that yields the following complete solution when the initial condition vanishes: 



(69) 



J to 



tr 



(70) 
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The one point two times correlation function then reads 

■>mm{t,t') / i \ ^7 



/mm{t,t') /. \ a7 

G-\T)i^-^j dr, (71) 

and after inverting Fourier we arrive at the real space expression 

{h{x,t)h{x,t')) = V(/i„(t)/i„(t'))cos2 , (72) 

where we have assumed no flux boundary conditions as in [l^, although the values of 
both the growth and auto-correlation exponents do not depend on the choice of boundary 
conditions. The propagator Gn{t) suggests the scaling variable Vn ~ nt^^"'^'^^/^ in Fourier 
space, that corresponds to the real space scaling variable u ~ as can be read 

directly from Eq. fl72l) . This suggests the definition of the effective dynamical exponent 
ZeS = C/(l ~ 7^)- If we express the correlation Eq. (17T|) for t = t' in terms of the scaling 
variable Vn (and we refer to it as C(f„) multiplied by a suitable power of t) and we introduce 
the "differential" 1 = An ~ t^'^+^^^^Av, we can cast the last expression in the integral form 

{hix, tf) - {h{x, t))' = ti-^/C+7<^WC-i) r c{v^) cos^ ( dv^, (73) 



where the series converges as a Riemann sum to the above integral when 

Dt«(4 + Z}to)^, (74) 
*o 

or equivalently t ^ tc ~ -^0°*^; for tf. being the time it takes the correlations reaching the 
substrate boundaries, assuming that the substrate initial size is very large. If 7 < 1/5, the 
whole substrate becomes correlated, yielding a finite tc-, for 7 > 1/5 the convergence of the 
Riemann sum to the integral is assured for all times, corresponding to the physical fact 
that the substrate never becomes correlated. In front of the integral we find a power of the 
temporal variable compatible with the growth exponent 

and the integral can be shown to be absolutely convergent as the integrand decays faster 
than exponentially for large values of the scaling variable f„. 

We are now in position to calculate the temporal auto-correlation 

,A _ {h{x,t)h{x,t'))o f min{t, t'} \ ^ 

^ - {h{x,mr{h{x,ml^' ^ [m^^{t,t'}) ' ^""^ 
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where A is the auto-correlation exponent and {■)o denotes the average with the zeroth 
mode contribution suppressed, as in fl73l) . The remaining ingredient is the correlation 
{h{x,t)h{x,t'))Q. Going back to Eg. (172]) we see that the Fourier space scaling variable now 
reads 



-7<5 



1/C 



1-,* J P^' 

If 7 < 1/6 the term max{t, t'}^""^^ is dominant and the factor in front of the convergent 
Riemann sum reads 

max{t, t'} (VC-i)7rf-rf/C inin{t, t'}, (78) 

after the time integration has been performed and in the limit max{t,t'} ^ min{t,t'}. 
In this same limit, but when 7 > 1/6, the term min{t, t'}^^'''^ becomes dominant and the 
prefactor reads 

max{t, t'}-"^^ mm{t, fji-^/fW/C, (79) 

The resulting temporal correlation adopts the form indicated in the right hand side of fITBl) . 
where 

^^h + d/C + ^d{l-6/0 if ^<l/6, ^^^^ 
[(3 + 'yd if 7 > 1/6, 

or alternatively 

A = /? + -, (81) 
where the A— dynamical exponent is defined as 

1/7 if 7 > 1/6. 



If we disregarded the effect of dilution we would find again Eq. fl8T|) . but this time 

= ZeS if 7 < 1/5, , ^ 

00 if 7 > 1/5. 

To further clarify the dynamics we now calculate the scaling form that the two points 
correlation function adopts for short spatial scales |a; — a;'| ^ t(i^<5')')/C jj^ the decorrelated 
regime. As dilution does not act on such a microscopic scale, the following results are 
independent of whether we contemplate dilution or not. In this case one has 

{h{x,t)h{x',t)) = 2^ '''^^m^l;^) ^°^( "lt^ 

(84) 
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where x = (xi, ■ ■ ■ , Xd) and n = (ni, ■ ■ ■ , n^), and we assume the rough interface inequahty 
( > d in order to assure the absolute convergence of this expression. By introducing the 
scahng variables Vi = riit^^'^"''^^'^ and Ui = Xit^'^^~^^^'^ for i = 1, - ■ ■ ,d and assuming statistical 
isotropy and homogeneity of the scaling form we find 



{h{x, t)h{x', t)) - {h{x, t)f = \x- x'\^-H^^^-^^T [\x - x'lt('^^-i)/^] 



(85) 



or in Lagrangian coordinates \y — y'\ = \x — x'lf^ 

{h{y,tMy',t)) - {h{y,t))' = |y - -^t^(^-^)^ 



\y-y 



t{l+7(C-5)}/C 

We see that this form is statistically self-affine with respect to the re-scaling y^by,t 
and h ^ b°'h, where the critical exponents are 



a 



C-d 



+ 



c 



(5 - C)7 



c 



(86) 
(87) 



Note that the scaling relation a = (3z holds, where the growth exponent (3 was calculated in 
Eq. fl75|) . The macroscopic decorrelation, which is observed for length scales of the order of 
the system size |x — x'| ~ Lq, is controlled by the effective dynamical exponent Zcs- When 

5 > C decorrelation might happen at microscopic length scales |x — x'| <^ t(i^^')')/C ^ell. 
Microscopic decorrelation happens in the limit 5 ^ C + I/7. For 5 < C + I/7 the interface 
is microscopically correlated and the critical exponents take on their finite values given in 
Eq. flHTj) . For 6 > C + I/7 the interface is microscopically uncorrelated and the critical 
exponents diverge a = z = 00, while the growth exponent is still finite and given by Eq. 
( l75l) (so one could say the scaling relation a = f3z still holds in some sense in the microscopic 
uncorrelated limit). With respect to the growth exponent we can say that (3 < 1/2 when 

6 < 7~^ + C) /5 ~^ 1/2 when 6 J'^ + Cy and P > 1/2 when 6 > 7~"'^ + C) so rapid roughening 
is a consequence of microscopic decorrelation. And now, by applying the developed theory 
to the IS equation, for which d = 2,( = 2,6 = 3 and assuming as in 15|] that 7 = 1, we find 
that it is exactly positioned at the threshold of microscopic decorrelation, this is, its critical 
exponents are a = z = 00 and /3 = 1/2. 

Note that the effective dynamical exponent Zes = (/{I — 1^) states the speed at which 
both correlation and decorrelation occur. The transition from correlation to decorrelation is 
triggered by the comparison among the indexes 7 and 6. The derivation order ( controls the 
speed at which both processes happen: a larger ( implies slower correlation/decorrelation 
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processes. Note also that rapid roughening might appear in exactly the same way in planar 
processes, just by allowing field or time dependence on the diffusion constant. This could 
be thought as somehow artificial in some planar situations, but as we have shown it appears 
naturally in the radial case, where such a dependence is a straightforward consequence 
of the lost of translation invariance, due to the existence of an absolute origin of space, 
characterized by a zero radius (and which in turn implies the existence of an absolute origin 
of time in the small noise approximation, as we have already seen). Such a naturalness can 



be seen in the derivation of the IS equation in [l5|, where it was found as a consequence of 
a simple variational principle. 

VI. POLYFRACTALITY 

We devote this section to showing that rapidly growing radial interfaces develop "polyfrac- 
tality" . We coin this term to denote a behavior characterized by a scale dependent fractal 
dimension taking place in a finite system and for long times. It is different from the concept 
of multifractality, which in this topic is usually associated to a nonlinear relation among the 
exponents characterizing the higher order height difference correlations 

In the classical case of static planar interfaces the fractal dimension is computed from 
the height difference correlation function 

{[h{x, t) - h{x\ t)fY'^ ~ |x - x'l^, (88) 

in the long time limit, i. e. after saturation have been achieved, where the Hurst exponent 
H = (C — c?)/2 for linear growth equations and the right hand side is time independent. The 
interface fractional dimension is calculated using the box counting method and is given by 
df = 1 + d — H. The general linear equation for stochastic growth on a growing domain was 
found in the last section to be 

dth = -D \V\^h -^h + ^Ff'-' + ^jy^\{x,t), (89) 

for which we will assume C < 6 < ( + 7^^. Its Fourier transformed version, for n > 1, is 

''^ -D ( ^1 -'-^h.+ ( t) '''' Ut). (90) 



dt \t Li t \t 
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For slow growth 7 < 1/5 diffusion dominates over dilution and one finds an expression 
compatible with that of the planar case 

{[h{x,t) - h{x',t)fY'^ ~ - x'\^'^-'^/\ (91) 

and so the Hurst exponent and interface fractal dimension are the same as in the planar case 
for fixed time. In the case of fast growth 7 > 1/5, for small spatial scales |x — x'| ^ t(i^'5'>')/C 
we recover again this result, while for large spatial scales |a; — x'l ^ t(i~^')')/C 

{Hx,t)-h{x\t)fY'' ^t^ (92) 

and so, for fixed time, H = and df = d + 1. This means that the interface becomes highly 
irregular and so dense that it fills the {d + 1)— dimensional space. This way decorrelation 
marks the onset of polyfractality, as specified by a scale dependent Hurst exponent, whose 
asymptotic values are 

Hi\x-xlt) = {^' (93) 
[0 if |x - x'\ > 

and the corresponding asymptotic values of the scale dependent fractal dimension 

, l + (3d-C)/2 if 
^ d + l if \x-x'\ 

Note that these results imply dynamic polyfractality as the scale separating the two regimes 
depends on time |x — a;'| ~ 

t{i-57)/C. also, the rough interface inequality C, > d implies the 
strict inequality 1 + {?>d — Q/2 < d + 1. This asymptotic behavior strongly suggests the 
self-similar form of both Hurst exponent and fractal dimension 

According to this the fractal dimension would be a dynamic fractal itself, invariant to the 
transformation x — > hx, t —* b^H, and df — > b"^df, for Zf = (/{I — ^7) = Zes and a/ = 
0. Note that all these results concerning polyfractality are independent of whether we 
contemplate dilution or not (because the height difference correlation function depends on 

n 

strictly local quantities PJ[), and so we could, in this particular calculation, substitute 
Eqs. ( l89l) and ( l90l) by their dilution-free counterparts and still get the same results. Note 
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also that at the very beginning of this section we have assumed the inequahty ( < S < C+7~^) 
which imphes that for rapid growth the interface is macroscopically but not microscopically 
uncorrelated. If 6 > C + 7~^ then the interface is microscopically uncorrelated and the fractal 
dimension becomes df = d + 1 independently of the scale from which we regard it, i. e., 
polyfractality is a genuine effect of macroscopic decorrelation, which disappears for strong 
damping causing microscopic decorrelation. 

Note that polyfractality does not appear in non-growing domain systems as for long times 
saturation is achieved and the fractal dimension becomes constant (assuming no multifrac- 
tality is present). Although the behavior of the height difference correlation function we 
found here is similar to the one present in classical unbounded systems, results concerning 
the fractal dimension cannot be immediately extrapolated. The fractal dimension can be 
computed in a bounded growing domain, for instance using the box counting method as 
we have done herein, by employing as the reference length L{t), the linear time dependent 
size of the system. Of course, in an unbounded static domain there is not such a reference 
length. 



VII. THE KARDAR-PARISI-ZHANG EQUATION 

One of the most important nonlinear models in the field of surface growth is the Kardar- 
Parisi-Zhang (KPZ) equation 25 1 



dth = uV'^h + \{\/hy + ^{x,t). (96) 

It is related to the biologically motivated Eden model, as this model, at least in a planar 
geometry, was numerically found to belong to the KPZ universality class As we will see, 
understanding the KPZ equation on a growing domain may shed some light on some of the 
properties of the classical version of this model. 
The KPZ equation on a growing domain reads 

d.h = V (^) V^/. + \ {^hf - + 7Ft-^ + '''' e(x, t). (97) 

Of course, if we just considered the dilation x — > {t/tQ)'^x we would find 

dth = V V'h + \ {Vhf + ^Ft^-' + e(x, t). (98) 
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As we have shown in the previous section, the dilution mechanism fixes the Family- Vicsek 
ansatz in the fast growth regime. In the radial Eden model case, assuming it belongs to 
the KPZ universality class, we would have z = 3/2 in d = 1 and 7=1. And so, one 
would naifly expect that the resulting interface is uncorrelated and we have to resort on 
dilution effects in order to fix the Family- Vicsek ansatz and get rid of memory effects. But 
here comes the paradoxical situation. There are two main symmetries associated with the 
(i-dimensional KPZ equation: the Hopf-Cole transformation which rnaps it onto the noisy 
diffusion equation 26|] and the related directed polymer problem |27|, l28|, and Gahlean 
invariance which have been traditionally related to the non-renormalization of the KPZ 

In the case of the 



vertex at an arbitrary order in the perturbation expansion [29 
no-dilution KPZ equation (jHSD both symmetries are still present. Indeed, this equation 
transforms under the Hopf-Cole transformation u = exp[A/i/(2z/)] to 



dtu 



27 



^2 7^^ ^ 1 ^ 



(99) 



2u 2u 

which is again a noisy diffusion equation and it can be explicitly solved in the deterministic 
limit e = 0. We find in this case 

|a;-y|2(l-27) 



u{x, t) 



(1 -27)^^/2 exp[FAr/(2z/)] 



[47rt^^(ti-27 -4-2^)]<i/2 
which corresponds to 



exp 



4t^^(ti-27_4-27) 



u{yM)dy, (100) 



2u 

h{x, t) = — In 
A 



;i-27)"'/^exp[FAt^/(2;/)] 
[47rt^^(ti-27 - tJ"'^)]^/2 



exp 



\x 



27) A 



dy 
(101) 



for given initial conditions u{x,tQ) and h{x,tQ). It is clear by regarding this formula that 
decorrelation at the deterministic level will happen for 7 > 1/2. It is still necessary to 
find out if at the stochastic level this threshold will be moved to 7 > 2/3. If we consider 
the dilution KPZ equation ( j97l) then transforming Hopf-Cole we would find the nonlinear 
equation 



dtu 



V 
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V\-^nln(n) + ^rt + A 
t ^ ' 2v 2v 



d-y/2 



^{x,t)u, 



(102) 



which may be thought of as a time dependent and spatially distributed version of the Gom- 
pertz differential equation 3l| . In this case it is not evident how to find an explicit solution 
at the deterministic level and what would be its decorrelation threshold. 
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Galilean invariance means that the transformation 

a; — > x — Aft, h —> h + vx, F ^ F — (103) 

where v is an arbitrary constant vector field, leaves the KPZ equation invariant. In case of 
no dilution this transformation can be replaced by 

A .9a,.i_o„, , . „ „ A 



x^x-- — h^h + vx, F ^ F ~ —vhl^t'-^\ (104) 

1 — 27 27 

which leaves invariant equation fl98l) . If we consider dilution, then it is not clear how to 
extend this transformation to leave equation (!97j) invariant. The main difficulty comes 
from the dilution term which yields a non- homogeneous contribution to the dynamics as 
a response to the rotation h ^ h + vx. So in summary we may talk of a certain sort of 
Galilean invariance which is obeyed by the no-dilution KPZ dynamics (!98l) and is lost when 
dilution is taken into account. If it were found that the dilution equation fl97j) obeys the 
traditional KPZ scaling (at least in some suitable limit), then that would mean the possible 
necessity for readdres sing the role that the symmetries of the KPZ equation have in fixing 



the universality class 32 



33, 



3^. 



There is still another fundamental symmetry of the KPZ equation, but this time it just 
manifests itself in one spatial dimension: the so called fiuctuation-dissipation theorem [l|. 
It basically says that for long times, when saturation has already being achieved, the non- 
linearity ceases to be operative and the resulting interface profile would be statistically 
indistinguishable from that created by the EW equation. For fast domain growth, we know 
from the linear theory that the interface never becomes correlated, and it operates, in this 

n 

sense, as if it were effectively in the short time regime for all times [16|. As a consequence, 
the fiuctuation-dissipation theorem is not expected to play any role in this case. Of course, 
this result would be independent of whether we contemplated dilution or not. 

In more general terms, it is known that the different symmetries of statistical mechanical 
models influence their scaling properties 35|, l36|]. It would be interesting to understand 
in complete generality the interplay among the symmetries of a physical model in a static 
domain and the asymmetric presence of dilution when we let this domain grow in time. 
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VIII. CENTER OF MASS FLUCTUATIONS 



Another property that has been studied in the context of radial growth, particularly in 
Eden clusters, is the center of mass fluctuations. It was found numerically that the Eden 
center of mass fluctuates according to the power law Cm ~ f^^^ in d = 1 + 1 {4], while in 
d = 2 + 1 there is a strong decrease in this exponent l37il . This reduced stochastic behavior 

n 

in higher dimensions was already predicted in [IJ] using radial growth equations, and we 
will further examine herein the compatibility among the equations and the Eden cluster 
dynamics. The center of mass fluctuations are characteristic not only of radial growth but 
also of planar situations. Let us recall the classical EW equation 

dth = D\/^h + ^{x,t), (105) 

defined on a one dimensional domain of linear size Lq and with no flux boundary conditions. 
It is straightforward to find that the center of mass ho{t) = Lq^ J^" h{x,t)dx is a Gaussian 
stochastic process defined by its two first moments 

{ho{t)) = 0, {ho{t)hois)) = min(t, s), (106) 

and so we have found that the center of mass performs Brownian motion, or equivalently 
we would say that its position is given by a Wiener process. Note that the fiuctuations 
amplitude decreases with the linear system size, suggesting that in the case of a growing 
domain our current law Cm = (^0)^''^ ~ ^^^"^ ^^^^ replaced by a different power law with 
a smaller exponent. It is easy to see that this result does not hold uniquely for the one 
dimensional EW equation; indeed, for any rf— dimensional growth equation with a conserved 



growth mechanism, be it linear as i\ie_ 
the Villain-Lai-Das Sarma equation 
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W or MuUins-Herring equations pL] or nonlinear as 



391] or its Monge- Ampere variation 23(], the center 



of mass performs Brownian motion characterized by the correlators 



(/io(t)) = 0, {h^{t)ho{s)) = -jmin(t,s), (107) 
as a consequence of the decoupling of the zeroth mode with respect to the surface fluctua- 



tions 



38 ] . Note that in the case of non-conserved growth dynamics this is not the case, as 



illustrated by the KPZ equation 

dth = vV^h + \{yhf + i{x,t). (108) 
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It is easy to see that in this case 

^ = ^1 i^^hfdx + Ut) > Ut), (109) 

where ^o(^) = '^J ^{x,t)dx and the equal sign is attained only for h = constant, an 
unstable configuration for KPZ dynamics. And so one expects stronger center of mass 
fluctuations in this case. Actually, the short time center of mass fluctuations can be easily 
calculated for any model which obeys the Family- Vicsek scaling, including the KPZ equation. 
Indeed, the Family- Vicsek scaling implies the following form of the height-height correlation 

{hix, t)h{x\ t)) = t'^ C , (110) 

which in the short time limit tends to 

{h{x, t)h{x', t)) ~ t'^+'^Z^ 6{x- x') , (111) 

leading to the result 

{ho{tf) ~ L-H^^+'^'\ (112) 

And so, within the Family- Vicsek scaling framework, the exponent characterizing the short 
time behavior of the center of mass fluctuations is /3 + d/{2z). 

As we have seen, the center of mass fluctuations are given by the zeroth mode. In the 
growing domain case it can be shown that the equation controlling the evolution of ho is 



16| 



^^^n, + ,Fr-^ + C-^T'''ut), (113) 



dt t \t ^ 

in case dilution is taken into account. In this case we find for long times the center of mass 
fluctuations 

Cl = {hom - {h,{t)f = JT^^f-'^ (114) 
and so Cm ~ t^^~'^'^'^/'^ . If we did not consider dilution we would find in the long time limit 

Cl=\ ^ln(t) if 7 = 1M (115) 

If we adapt result flll2p to the present setting we find 

{h{x, t)h{x\ t)) ~ 5 [i7(a; _ a;')] = ^2p+d/z^^d ^r^ _ ^^^g) 

25 



For linear systems the equality 2[3+d/z = 1 holds, and so this last equation agrees with flll4p 
but not with flllSp . This is a consequence of the violation of the Family- Vicsek scaling in the 



absence of dilution 
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18|. In the case of the (1 + 1)— dimensional Eden model c? = 7 = 1, 



and if it belonged to the KPZ universality class the center of mass would fluctuate according 
to the law Cm ~ t^^^. This of course does not agree with the measured behavior Cm ~ t"^^^- 
This exponent could be recovered by introducing an ad hoc instability mechanism, such as 
for instance considering a growth equation whose zeroth moment obeyed 

57 /+ \ d.'y/2 



dho 
~dt 



D(^jyho + ^Ft^-'+(^jy Ut). (117) 



The desired exponent is obtained for 6 = 1 and Dto = 2/5, but however this result is uniform 
on the spatial dimension and so can not predict the (2 + 1)— dimensional behavior js^]- 
Additionally this instability mechanism seems to be not enough justifled and too non-generic 
to be a good explanation of the observed phenomenology. Everything points to the fact that 
the center of mass fluctuations of the Eden model result from a strong violation of the Family- 
Vicsek ansatz. As we may see from equation flllSp . this sort of violations imply stronger 
center of mass fluctuations. This point will be further discussed in the next section. 

In summary we can say that the result Cm ~ t"^^^ suggests a strong violation of the 
Family- Vicsek scaling by the surface fluctuations of the (1 + 1)— dimensional Eden model. 
Although the linear law Cm ~ ^(^^t'^)/^ does not reproduce quantitatively the results, we 
still expect from it a qualitative description of the dynamics, as the strong decrease of this 
exponent was already reported in (2 + 1)— dimensions. According to the linear law, the center 
of mass fluctuations should decrease for increasing growth velocity and spatial dimension. 
Note also that the nonlinearity seems to be a necessary ingredient; the linearization of the 



KPZ equation proposed in 



ll| reads in Fourier space 



±{hl) ^ -A\n\"^kl) + (118) 

for some constants A and B and in case of a non-growing domain. This equation supports 
unbounded fluctuations as revealed by the divergent diffusion in the limit n —>■ 0, and so 
this does not constitute a good model for predicting the center of mass fluctuations. 
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IX. APPLICATIONS TO THE EDEN MODEL 



In statistical mechanics it has been customary to classify the behavior of discrete models 
within universality classes defined by continuum field theories. Non-equilibrium growth the- 
ories have been by no means an exception to this rule 0, |4l|. In this sense one would be 
interested in finding the universality class the Eden model belongs to. According to the sim- 
ulations performed in the planar geometry the Eden model belongs to the KPZ universality 
class [rl. This agrees with the measured exponent /5 = 1/3 in radial systems ^j. However, 
as we have already seen, there are at least two possible universality classes associated to the 
KPZ equation in radial systems: dilution-KPZ and dilation-KPZ. The first one is character- 
ized by a behavior more akin to that of planar systems, and the second one by memory effects 
which imply the departure from the Family- Vicsek ansatz. According to the measurement 
of the autocorrelation exponent of the Eden model in ll|] that yielded A = 1/3, the Eden 
model would be in the dilation-KPZ universality class (one would expect A = 4/3 for dilution 
dynamics according to the theory developed herein and in [l6|). This fact admits a simple 
explanation. In the Eden model, cells are aggregated to the colony peripheral in such a way 
that the positions of already present cells are not modified. Consequently, as the system 
grows, no dilution is redistributing its constituents. So the rigidity of the Eden model may 
well be at the origin of the memory effects present at its interface jll|, which presumably 
place it in the dilation-KPZ universality class. But to be sure one would still need, of course, 
to verify that this implies no contradiction with the center of mass fluctuations as discussed 
in the last section. 

As we have already mentioned, the Eden model may be thought of as an idealization of a 
developing cell colony. Of course, as it was completely clear from the very beginning 
there are multiple factors of biological, chemical and even physical nature that are not 
captured by this model. Apart from them, one could be interested in improving the model 
in pure statistical mechanical terms. To this end one may look for inspiration in real cell 
colonies. The structure of a rapidly developing cell colony would be dominated by dilution 
effects, originated in the birth of new cells which volume causes the displacement of the 
existent cells. This feature is not captured by any sort of Eden model (diverse proliferation 
rules, on/off lattice,...) and is fundamental in preserving the Family- Vicsek scaling, as we 
have already seen. So it seems quite reasonable to modify the Eden model in order to 
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remove its rigidity, allowing bulk cells proliferation and the displacement of the existent 
cells, both at the bulk and interface, by the newborn cells. This would not be interesting 
just in modelling terms, but also for introducing dilution in the model and consequently 
shifting its universality class. 

X. CONCLUSIONS AND OUTLOOK 

In this work we have investigated the role of dilution and decorrelation on radial growth. 
Dilution drives matter redistribution along the growing interface: as the surface becomes 
larger the already deposited matter occupies a smaller fraction of interface, which is being 
simultaneously complemented with incoming matter, the actual driving force of domain 
growth in radial systems. Dilution is important for any rate of domain growth, as it keeps 
the interfacial density constant, but specially for rapidly growing domains, for which the 
diffusion mechanism becomes irrelevant and dilution becomes the sole responsible for the 
propagation of correlations on the macroscopic scale. The importance of dilution is such that 
in its absence, which takes place in the alternative dilation dynamics, strong memory effects 
arise. These include an enhanced stochasticity, which separates the behavior of the large 
spatial scale limit of the two-points correlation function from that dictated by the Family- 
Vicsek ansatz, and the appearance of non-universal critical exponents in the marginally 
rough regime, characterized by the equality C — d. As have seen, both universality and the 
Family- Vicsek structure of the correlation function are recovered by virtue of dilution. 

As dilution propagates correlations at the same speed at which the interface grows a 
global correlation becomes impossible for fast domain growth. This leads to decorrelation, 
or in other words, to a whitening of the interfacial profile in the sense that distant points 
become uncorrelated. Decorrelation might be macroscopic, which is evident only if we 
regard the dynamics from a spatial scale of the same order of magnitude of the system size, 
or microscopic, in which case it is apparent for much smaller length scales. Microscopic 
decorrelation supports rapid roughening, i. e., growing regimes characterized by (3 > 1/2. 
These appear naturally in the context of radial growth, for instance by considering the IS 
equation, which results from a geometric variational principle and for which = d = 2 and 
(5 = 3, and thus it shows rapid roughening for all 7 > 1. A consequence of macroscopic 
decorrelation is the advent of a scale dependent interfacial fractal dimension, which rends 
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the surface polyfractal and we have conjectured it to be self-similar. 

There are several theoretical problems that can be straightforwardly analyzed with the 
techniques introduced here. We have for instance considered radial interfaces whose mean 
radius grows as a power law of time (r) ~ t"'. This result has been obtained by means of a 
linear mechanism in which an explicit power law dependence on time has been considered, 
see Eq. (fTOj) . This linear mechanism can be substituted by a nonlinear one in which time 
does not appear explicitly 

a,. = -,F'/V-'/^+-^^^«9.t), (119) 

which yields at the deterministic order R = FV again, but it is the source at the first stochas- 
tic order of a term (reminiscent of dilution) which may be either stabilizing or destabilizing 
depending on the value of 7 

dtp=^^^p + — - — (120) 

for small values of 7 the previous sections results are recovered, while for large values of 7 
memory effects and enhanced (power law) stochasticity appear (which are standard effects 
of instability as we have already seen), with the threshold value of 7 depending of whether 
we introduce dilution or not (in this concrete example dilution completely erases instability). 
Also, this instability mechanism, contrary to the ones studied herein and in which rend 
the zeroth mode unstable and the / = 1 ones marginal, is able to destabilize all modes. 
Different nonlinearities which might destabilize a fixed number of modes lying before some 
given Z* G N can be easily devised too (basically by introducing terms of the form — r~™ 
for some suitable m G N in the corresponding equation of motion) and can even be cast on 
some geometric variational formulation as the cases considered in [l^. Of course, deciding 
which model is the good one must rely on numerical or experimental evidence based on the 
study of specific models or systems of interest. 

As mentioned in the introduction, part of the motivation for studying radial growth 
models such as the Eden or different ones lies in the possible similarity of these with some 
forms of biological development, such as for instance cell colonies. The results of our study 
can be translated into this context to obtain some simple conclusions, provided the modelling 
assumptions make sense for some biological system. The structure of a rapidly developing 
cell colony would be dominated by dilution effects, originated in the birth of new cells which 
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volume causes the displacement of the existent cells. If the rate of growth is large enough 
this motion will dominate over any possible random dispersal of the individual cells. It is 
remarkable that such a consequence simply appears by considering domain growth, while 
it is not necessary to introduce corrections coming from the finite size of the constituents. 
This is the dilution dominated situation we have formalized by means of the (decorrelation) 
inequality 7 > 1/C (assuming in this case S = (). If we were to introduce some control 
protocol in order to keep the consequences of bacterial development to a minimum we would 
need to eliminate colony constituents (possibly randomly selected) at a high enough rate 
so the effective growth velocity were one that reversed the decorrelation inequality. For the 
one dimensional Eden model, accepting it belongs to the KPZ universality class, one finds 
7 = 1 and z = 3/2. If z played the same role for the nonlinear KPZ equation as ( for the 
linear equations considered herein (as it is reasonable to expect), the Eden model would be 
in the uncorrelated regime. In order to control it we would need to eliminate its cells at rate 
such that the effective growth rate obeyed 7 < 2/3. For the two dimensional Eden model, 
if its behavior were still analogous to that of the KPZ equation, we would find z > 3/2 and 
thus a greater difficulty for control. Note that for the particular growth rules of the Eden 
model one would need to eliminate peripheral cells in order to control the system. This 
would not be so in the case of an actual bacterial colony, for which bulk cells are still able to 
reproduce, and so cell elimination could be performed randomly across the whole colony. Of 
course, these conclusions are speculative as long as radial growth equations are not proved 
to reasonably model some biological system. 

In more general terms, we have found that the surface fluctuations of the Eden model 
presumably strongly violate the Family- Vicsek scaling. We have identified the absence of di- 
lution in this model as the reason underlying such a violation. In this sense, this model would 
not be able to describe growing cell colonies, precisely because it assumes a spurious rigidity 
of bulk cells. On the other hand, it would be better suited to describe the radial growth 
of crystalline structures 42]. We have also found that reparametrization invariance [2^ 
implicitly implies dilation dynamics. Our results call for an extension of the generalization 
of Langevin dynamics to arbitrary geometries in order to capture both dilution and dilation 
scenarios, and the associated bifurcation of universality classes. This same remark would 
affect as well equilibrium systems, but in this case of course the domain evolution will drive 
them out of equilibrium, unless growth is quasistatic j^sj]. 
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APPENDIX A: HIGHER ORDER PERTURBATION EXPANSION 

As we have mentioned in Sec. [TTl the first order correction in the small noise expansion 
is a Gaussian stochastic process. We will try to go beyond this order in this appendix, and 
we will show the difficulties that arise in trying so. We focus again on the radial random 
deposition equation (fTOj) and assume the solution form 

r0,t) = R{t) + V~ep{lt) + ep^ilt), (Al) 

where the noise intensity e will be used as the small parameter |21j] . Substituting this solution 
form into Eq. ffTOj) we obtain the equations hierarchy 

dtR = F^r-\ (A2) 
dtp, = "^^t'^^ , (A3) 

where ^ = ^/t t] and both r] and ^ are now zero mean quasiwhite Gaussian processes whose 
correlations are given by 

r]{lt)i^{lt)) = c{e-6')5{t-t'), U(e,me,t)) = eC{e-e')6{t-t'), (A5) 



where C(-) is some regular function approximating the Dirac delta; the necessity for the 
quasiwhite assumption will we clear in few lines. These equations have been derived assum- 
ing ^/e <^ Ff^, and we will further assume a zero value for both initial perturbations as in 
Sec. [TTl The solution to the first two was characterized in Sec. [TTl where the approximating 
function C(-) was substituted by the Dirac delta. Here i? is a deterministic function and pi 
is a zero mean Gaussian stochastic process that is completely determined by its correlation 
function. The stochastic function p2 is a zero mean process too, but it is not Gaussian this 
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time, and its correlation (which no longer completely determines the process) is given by 



P2{e,t)p2{e',s] 



X 



mm 



in{t, s})2-2T-2T<i _ f-2j-2jd ^^-^^^^ sjy-2j-jd _ ^1-27-7^ 



2 - 27 - 27d " 1 - 27 - 7d 

if 7(i 7^ 1, 7(1 + d)j^l, and 7(2 + d) 1. li -fd = 1 we find 



/A2 



c{e - 9') 



; (A6) 



if 7(1 + d) = 1 then 
[p2{9,t)p2{9',s) 



to - [min{t, s}] 



-27 



1 + 27ln 



mm 



0,4 



\ C{9- 9'f 
(A7) 



4_p2+2d^ 

and if 7(2 + 0?) = 1 we get 
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min{t,g} \ ^ 
to ) 1 



+ ^([minO,s}]-^-io'') 
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P2(^"',t)p2(^"'',s)) = 



g^2+2d^ 



27 



tn^ln 



min{i, s} 
To 



C{9 - 9'f 
J{9)J{9') 



(A8) 



■• (A9) 



The long time behavior of the correlations, given by the condition t,s ^ to, is specified by 
the following two-times and one-time functions 

\ r/2 ^2-27-27d C(0 — f)'\2 

when 7((i -|- 1) < 1, and if 7(0? + 1) = 1 then 



P2{9,t)p2{9',s) 



-—7-— 3— in (mmjt, s\) =; —. 

4F2+2'^7 ^ ^ ' J{0)J{9') 



P2{9,t)p2{9',t) 

and finally, when 7(rf -|- 1) > 1, we find 
p^(lt)p.^(S',s)' 



4^2+2^7 J{9)J{9') ' 



(A12) 
(A13) 



2-27-27^ 




A 2 



(A14) 



8F2+2rf 1 - (3 + 2d)7 + (2 + 3(i + ^2)72 J{9)J{9') ' 

a correlation function that vanishes in the limit ^ 00. Now it is clear why we needed 
the quasiwhite approximation: for a regular function C(-) the expression C(-)^ makes sense. 
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contrary to what happens if we substitute it by the Dirac delta to get S{-)'^. This is the first 
indication of the failure of the higher order perturbation theory. 

We now examine the effect that dilution has on the random function p2, which in this 
case obeys the equation 

_ d (rf + i)w p,{e,mlt) 

In this case the long time correlation function reads 

pMS;s)) ^ ^) its)-' ■n.nt, (A16) 



if 7 < 1, 



if 7 



if 7 > 1. The one time correlation function is then 

if 7 < 1, 
if 7 = 1, 

'^^^''^'^^^''V-8F2+-(7C^ + l)(7-l)'° J(t-./)J(t-.f)' ^^''^ 

if 7 > 1, where £ - f = t^(^- C*!^"- ^'^) = ^'), and we have assumed that 

the approximating function C(-) has the same homogeneity as the Dirac delta. Although 
it is evident that dilution carries out a measurable action, particularly erasing part of the 
memory effects, the result is far from satisfactory. In all cases the prefactor deviates from the 



expected random deposition form 



44| . the unexpected critical value 7=1 has appeared, 



and for 7 > 1 memory effects are present as signaled by the logarithm and the to dependence 
respectively; and the situation is further complicated by the presence of the factor C(-)^ 
which becomes singular in the white noise limit. All of these elements suggest the failure of 
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the small noise expansion beyond the first order. Classical results suggest the possibility of 
constructing a systematic approach to the solution of some nonlinear stochastic differential 
equations by continuing the small noise expansion to higher orders 2l| . Our present results 
suggest the failure of this sort of expansions beyond the Gaussian (which turns out to be 
the first) order in very much the same way as the Kramer s-Moyal expansion of the master 



equation 



45l | and the Chapman- Enskog expansion of the Boltzmann equation [46| fail beyond 



the Fokker-Planck and Navier-Stokes orders respectively. 
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